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A magnetic flux line lattice or a vortex lattice, where each flux line contains a quantum 
of flux (^Q=hc/2e) and interacts repulsively, is expected to occur in the entire mixed 
state of a type-II superconductor, within a mean field description . Inclusion of thermal 
fluctuations predicts melting of this vortex lattice into a vortex liquid state, for both 
dilute and concentrated arrays of lines, yielding an unusual reentrant phase 
boundary^^ in the magnetic field-temperature (H,T) plane, i.e., at a fixed T, the melting 
phase boundary is encountered twice. Quenched pinning disorder adds more variety in 
the form of novel disordered (glassy) phases^ \ and makes the vortex lattice an ideal 
system to study the competition and interplay between interaction and disorder. Our 
experimental studies of magnetic response in the mixed state of the superconductor 2H- 
NbSe2 (in samples of progressively increasing pinning) show the appearance and 
evolution of novel behavior such as history dependence, reminiscent of spin glasses^ 
These results elucidate how pinning alters the phase boundary separating the 
"ordered" and "disordered" states of the flux lines and causes a novel disorder- 
induced transition. A generic magnetic phase diagram of a type II superconductor, 
that includes effects of thermal and quenched disorder, has been constructed. 



All the data are obtained using the standard magnetic ac susceptibility studies involving the 
response to a small amplitude ac magnetic field (hac) superimposed on a background dc bias 
field (Hdc) maintained parallel to the c-axis of two layered hexagonal Niobium diselenide 
(2H-NbSe2). The magnetic induction (B) determines the lattice constant 

a(,(=(24>(,/BV3)^'^ »5xlO^AV7fi(g«M*5))of the triangular flux line lattice (FLL)'. 
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Fig. 1 shows the typical variation of the real part of the ac susceptibility (% ) with temperature 

(T) : it shows perfect diamagnetism at low T -1) and rises to the normal state value at 
Tc(0) in zero-field. For a dc bias field, % can be approximated, in the Critical State modeF, 
by : x' ~ (-1 + ohjh), where a is a geometry dependent factor and Jc is the critical current 
density at a given (H,T). With increasing T, monotonically increases due to an expected 
decrease in Jc. As shown in the curve for H = 4 kOe in Fig. 1, this usual behavior in %' is 
interrupted at the onset of an anomalous enhancement of diamagnetic screening response ; 
reaches a sharp minimum at Tp ( which can be sharper than even the zero-field 
superconducting transition, as in Fig.l) , above which it recovers rapidly towards the normal 
state value at Tc(H). This nonmonotonicity in % is a consequence of a nonmonotonicity in J^. 
The minimum in %' thus corresponds to a peak in Jc , the ubiquitous peak effect^ (PE). In 
weakly pinned systems, Jc is given' by the pinning force ( JcB ) equation : JcB = 

2 1/2 

(np<fp >/Vc) , where Up is the density of pins, fp is the elementary pinning interaction, 
proportional to the condensation energy and Vc is the correlation volume ( of a Larkin"^ 
domain within which flux lines retain their order) directly related to the elastic moduli of the 
FLL. The PE results' from a rapid collapse in Vc of the Larkin domain ( more rapid than the 
decrease in <fp^>, which causes the usual monotonic decrease in screening response with 
increasing temperature or field) so as to produce an overall enhancement in Jc and hence in 
the diamagnetic response. This unusual collapse of Vc, i.e., the loss of order of the FLL, is 
triggered by an anomalously rapid softening of the FLL. Although precise quantitative 
explanations are stiU lacking, all available scenarios'" " intimately relate PE to the melting of 
the FLL. Indeed, the PE boundary in a few weak-pinning samples of 2H-NbSe2 were found'^ 
to mimic closely the theoretically predicted^ reentrant melting phase boundary of pinning free 
FLL. In what follows, we show in samples of varying quenched disorder (i.e., pinning 
centers) the evolution of the loss of order of FLL near the incipient melting transition. We 
focus on the temperature regime in between the onset of the peak effect and the 
reappearance, above Tp, of a conventional approach of to the normal state value ( hereafter 
referred to as the peak effect regime). Our results provide new insight on the details of the 
process of loss of order in the FLL and the appearance of phenomena such as history 
dependence ( and metastability ) and new pinning induced phase fransitions. 

I. Disorder, Metastability and History dependence at moderate fields : 
Three different single crystals (samples A, B and C, made from starting materials of 
progressively decreasing purity), spanning two orders of magnitude differences in critical 
current densities (Jc values in A to C range from few A/ cm^ to few hundred A/ cm^ at 
T=4.2K and in H ~ 10 kOe) were studied. In each of these crystals, the fiux lattice is well 
formed : the transverse correlation length Rc' is much larger' than the values of the lattice 
constant ao at moderate to high fields ( H > 1 kOe ). 

Fig.2 shows the results of a history-dependence experiment. Two different sample histories, 
commonly associated with disordered magnets such as spin-glasses, were investigated. In 
one, the sample is first heated to the normal state and cooled in the presence of a field (field- 
cooled, i.e., FC). In the second, it is cooled in zero-field and then the field is applied (zero- 
field cooled, i.e., ZFC). In both cases x'data were collected during the subsequent warm-up. 

Fig.2(a) shows the temperature dependence of %' for two values of the field, 15 kOe and 2 
kOe, in the cleanest sample A ( made from ultrapure starting materials, notably, nearly Ta- 
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free Nb). There is no measurable history dependence; the ZFC and FC curves are the same, 
i.e., the flux lines system reaches the ordered state regardless of the path chosen. In Fig. 2(b), 
the data from a typically clean crystal (sample B made from 99.9% pure starting materials) 
are shown for three values of the magnetic field. For H =15 kOe, the same field value as in 
Fig.2(a), history dependence appears : the FC state is slightly more diamagnetic at the lowest 
temperature. The onset of the peak effect regime occurs discontinuously through a 
downward jump in x at Tpi. The first peak in % at Tp is followed by the appearance of a new 
peak (seen as a shoulder on the rapidly declining background response in %'). The two 
curves, FC and ZFC, merge into one near this shoulder. For a lower field value of 5 kOe ( 
see inset in Fig. 2(b) ), both the features, the jump in at the onset of peak regime and the 
history dependence in it, disappear as in Fig. 2(a), although the two peaks spanning the peak 
regime are stiU identifiable ( cf. arrows in the inset of Fig. 2(b) ). At a further lower field of 2 
kOe, only one peak remains and the % curve now ( see inset in Fig. 2(b)) resembles the two 
curves in Fig.2(a). As stated above, at a given (H,T) value, the sample B has a larger Jc than 
that in the sample A. The appearance of structure and history dependence in ( in sample B ) 
in the peak effect regime thus appears to be a signature of enhancement of disorder / pinning. 
These results also suggest that in a given sample ( here, crystal B ), one produces an 
effectively stronger pinning, upon increasing H, as has been postulated earlier elsewhere^'^^ 

In Fig.2(c), we show the data for our most strongly pinned FLL (in sample C, made from 
commercial grade starting material). In this sample, even for H=5 kOe, a large history 
dependence appears. ( Note the similarity between the results shown in Fig. 2(c) and those 
obtained near the spin glass transition temperature (Tg) in a canonical spin glass system, like, 
CuMn , in which ZFC and FC magnetization curves merge at Tg/".) The FC state in Fig. 2(c) 
is more diamagnetic, i.e., more strongly pinned than the ZFC. This is equivalent to the 
transport measurements"* on another sample from the same batch of crystals of 2H-NbSe2, 
which showed a far greater critical current density for the FC branch and which could be 
annealed ( by a driving current ) into a more ordered and more weakly pinned ( smaller Jc ) 
state, Uke, the ZFC branch. The FC state in Fig. 2(c) results from the supercooling of the 
disordered vortex liquid state from higher temperatures; it nearly retains frozen-in liquid 
correlations and thus shows minimal signatures of any further loss of order for FLL on 
warming up across the peak regime. For the sample C, at lower fields ( < 2 kOe, data not 
shown here ), where pinning is effectively reduced, the situation is similar to that seen in 
sample B at 15 kOe. The history dependence and metastabUity are thus unambiguously 
related to increased pinning. 

The ZFC branch in sample C shows a very broad peak effect regime ( see Fig. 2(c) ). 
However, unlike results in sample B ( cf. Fig. 2(b) ), a clear multi-peak structure is 
unidentifiable in sample C. Instead, the onset of peak effect regime is through a large jump in 
% at Tpi, and ending also with a large jump at Tp, above which much of the history 
dependence in disappears. The sizes of the two jumps are amplitude dependent; they 
become larger with lower ac amplitude, implying that they are of an equilibrium origin, i.e., 
they mark two (discontinuous) phase transitions. Thus, with increased effective pinning, the 
PE feature of the flux lattice system evolves from a single history- independent and sharp peak 
( in sample A ) to a strongly history dependent broad peak effect region ( in sample C ), 
which is spanned by two discontinuous transitions for the more ordered ZFC branch, but is 
nearly featureless, without a conspicuous PE, for the more disordered FC branch. 
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We propose that Tpi marks a transition from a nearly defect-free lattice'^ such as a Bragg 
glass"', to a highly defective plastically deformed lattice"'''^''^ full of topological defects 
(dislocations) analogous to a "vortex glass" phase proposed by Gingras and Huse'^ This 
identification is supported by two other sets of experiments'^'^" on the same 2H-NbSe2 
system. In one'^, the onset of the peak effect was seen to be accompanied by plastic flow 
(spatially inhomogeneous motion) of a highly defective lattice. In the other case, noise 
measurements^" showed that the aforementioned defective moving state is closely related to a 
defective pinned state, implying that the plastic flow occurs through channels defined by 
dislocations that persist in the pinned state. 

Most of the history dependence in disappears (near) above the jump seen at Tp, i.e., the 
flux lines system loses order in equilibrium above Tp ( thus it could be well described by a 
pinned liquid^'^' or an amorphous lattice, implying the occurrence of FLL melting ). The 
correlation length Rc of flux lattice presumably becomes"* of the order of ao and remains nearly 
unchanged at higher temperatures. The collapse of the residual pinning at even higher 
temperature ( T > Tp ) is then caused by the rapid decrease of the pinning interaction fp. The 
inset of Fig. 2(c), comprising plots of Tpi and Tp vs H in sample C, further illustrates the role 
of increased effective pinning at higher fields through the broadening of the peak regime 
between Tpi and Tp and the concomitant shrinkage of the ordered state below Tpi. Tpi(H) and 
Tp(H) lines may indeed meet at a multicritical point at a low field value. 

Above Tp, the vortex lattice remains pinned upto at a higher temperature Tm, identified in 
the inset of Fig.l for sample A. ( Similar data in samples B and C, not highlighted here for 
brevity. ) Above Ti^r, % becomes positive, due to the differential paramagnetic effect^^ (DPE) : 
dM/dH > 0, for the reversible part of the M-H curve as diamagnetism decreases with increase 
in H (dM/dH ~ K ^ near Tc(H); for NbSe2 K ~ 10). Tm occurs near the upper edge of the 
peak effect regime and is frequency dependent, implying that depending on residual pinning, 
there is a further (dynamic) crossover from a pinned liquid state to an unpinned vortex liquid 
state as the viscoelastic relaxation time ( plastic time) decreases with increasing T^'^\ 

II . Effects of disorder at low fields : 

Now we concentrate on the low field regime. We recall that in this regime a reentrant 
behavior in Tp vs H was observed in the 2H-NbSe2 system'^, consistent with theory^ and 
simulations^\ This regime has been sparingly studied experimentally and little is known about 
its characteristics and the effect of pinning on it. Fig. 3 shows the variation of x' with T at low 
fields for the three 2H-NbSe2 samples under study. For all the samples, the zero- field y! 
curves have also been plotted for easy comparison. In contrast with the high field data - 
where increasing field leads to an effectively larger pinning - in the low field regime the 
apparent dominance of pinning with decreasing field is caused by the weakening of the inter- 
vortex interaction as the lattice constant ao approaches the range of interaction, which is 
comparable to the penetration depth A, . At low fields ( for H < 1 kOe ; ao > 1.5 x 10^ A°, 
indeed comparable to A, ( // c ) of 2.3 x 10^ in 2H-NbSe2'), one does not see the features 
observed above, e.g., neither a double peak structure nor jumps spanning the peak regime, 
etc., are present. Instead, one finds the general trends described below and summarized in the 
inset of Fig. 4 : 

(1) The peak regime broadens (i) with increasing pinning ( cf. Figs. 3(a) to 3(c) for the three 
samples ) and (ii) with decreasing field ( for a given sample ), as indicated by the error bars in 



4 



TIFR/CM/97/104(I) 



the inset of Fig. 4. This can be understood to result from an enhanced inhomogeneity in the 
sizes of Larkin'^ domains within a given sample at a given field and, thence, to a larger spread 
in the melting temperature of the entire flux line system. This would imply that for a highly 
disordered soUd, the melting phenomenon is a diffuse transition. 

(2) For the upper branch of the Tp curve, i.e., FLL melting curve, the Tp(H) boundary shifts 
down to lower T for more strongly pinned samples ( see, inset of Fig.4 ), showing that a 
disordered soUd ( which is correlated over a smaller volume ) requires less thermal 
fluctuations to melt it^^ 

(3) The reentrant lower branch of the peak effect (melting) line is seen clearly only for sample 
B and not so for the other two samples (see inset of Fig. 4). The Tp(H) data at low fields in 

the cleanest sample A shows a pronounced departure away from the upper melting curve 
well-described by the standard Lindemann melting relation^ : B^^ (ciVGi)H(;2(0) (l-T/TJ^ 
with |3™= 5.6, Gi = 3 X 10", He2(0) = 4.6 T ^ '^ and a Lindemann number Ci ~ 0.15. The 
observed behavior is in accord with a recent' theoretical scenario which shows a similar rapid 
drop of the melting curve, away from the Hc2(T) phase boundary, at low fields (i.e., above the 
nose or the turnaround of the melting boundary). The reentrant phase boundary may thus be 
at much lower fields in the very clean crystal A, as proposed theoretically^ and outside the 
precision of the present experiment This implies that, with increasing quenched disorder, 
the ordered phase is destabilized from above, on the upper curve, by thermal fluctuations 
and from below on the lower curve, by weaker interaction, thereby shrinking the ordered 
vortex solid phase region on both branches of the melting curve in the (H,T) phase diagram. 

(4) Finally, the disappearance of the peak effect altogether at lower fields is consistent with 
the scenario that a transition occurs between an entangled liquid phase of vortices to a 
pinning dominated disentangled liquid^ or a glassy phase at the Nelson - Le DoussaF'* line. 
This hne marks the locus ( in H ) where the entanglement length equals the pinning length' 

1/2 

and scales with (jc) ; this accounts well for the increased field values for the disappearance 
of PE peak in the more impure crystal. 

Taking aU these results together we construct , in the main panel of Fig.4, a generic phase 
diagram of a type II superconductor in the presence of both quenched disorder and thermal 
fluctuations. The nearly defect-free ( elasticaUy deformed ) vortex solid undergoes loss of 
order in steps. At the Hpi hne, an elastic-to-plastic transition leads to a plastically deformed 
solid which melts into a highly viscous pinned liquid at Hp. Depending on pinning, both the 
transitions ( at Tpi and at Tp ) can be discontinuous or continuous (within experimental error). 
At a higher field, Hm, a (presumably dynamic) crossover occurs into an unpinned liquid state. 
Hc2 marks the final crossover into the normal state. At low fields, Hpi and Hp approach each 
other ( see Fig. 4 ) and a single transition from the vortex soUd to a liquid below a 
multicritical point cannot be ruled out. At still lower H, and particularly near the nose of the 
peak effect phase boundary, the melting transition broadens significantly, eventually 
becoming undetectable as a peak effect. This experimentally ascertained phase diagram can 
be compared with theoretical expectations ( e.g., Fig. 25 of Ref. 5 ) for the mixed state. 

It is important to note that, as shown above, the location and nature of these phase 
transitions and crossovers are explicitly dependent on effective pinning which produces a 
richness of behavior such as metastability, history dependence and disorder induced 
transitions not present in the cleanest system. Our results demonstrate that these effects are 
experimentally accessible and call for more theoretical as well as experimental work. 
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Figure Captions : 

Figure 1 :Real part (x') of the ac susceptibility response [hac = 1.0 Oe (r.m.s) and f = 21 1 Hz] 
in our purest (sample A) crystal of 2H-NbSe2 in ( nominally ) zero bias field and in Hdc of 4 
kOe, respectively. The curve in Hdc = 4kOe shows a negative peak at Tp (i.e., peak effect 
(PE)), prior to Tc(H), whose sharpness exceeds that of the zero field superconducting 
transition. The inset reveals a small paramagnetic peak due to differential paramagnetic effect 
(DPE) in between the irreversibility temperature Tirr(H) and Tc(H). 

Figure 2 : %'(T) curves for field cooled (FC) and zero field cooled (ZFC) states at selected 

fields in three crystals of 2H-NbSe2. ZFC data points are marked in red whereas continuous 
lines (in blue) pass through closely spaced FC data points (not marked). In cleanest sample A 
(panel (a)), FC and ZFC data points overlap in aU fields, whereas in samples B and C (panels 
(b) and (c)), the differences between FC and ZFC data can be seen in H = 15 kOe and 5 kOe, 
respectively. Tpi and Tp, respectively, mark the temperatures of the onset of anomalous 
change and the maximum of diamagnetic screening response in in each of the figures. The 
insets in Fig. 2(a) and Fig. 2(b) show FC and ZFC data recorded at H = 2 kOe in sample A 
and at H = 5 kOe and 2 kOe in sample B, respectively ( see text for description). The inset in 
Fig.2 (c) shows variation of Tpi and Tp with H in sample C. 

Figure 3 : x'(T) curves at different H ( < 1 kOe ) in crystals A, B and C of 2H-NbSe2 . The 
zero field curve is sharpest in A and broadest in C, thereby showing the progressive decline 
in the purity. Also note, the peak regime broadens and Tp(H) (for a fixed H of 300 Oe ) shifts 
to progressively lower reduced temperature value on going from A to C. In sample A, the PE 
peak remains very sharp even at H = 30 Oe, whereas in sample B, it broadens rapidly below 
300 Oe. In sample C, a distinct PE peak is not identifiable below 300 Oe. 

Figure 4 : A schematic plot of magnetic phase diagram in a weakly pinned type - II 
superconductor ( here, for example, crystal B ). Vortex ( mixed ) state is sandwiched 
between the Meissner state ( below lower critical field Hd line) and the normal metaUic state 
( above upper critical field Hc2 ). Different phases of vortex lines have been identified on the 
basis of experiments in 2H-NbSe2 system. Hpi, Hn,eit, Hi^r and lines correspond to Tpi(H), 
Tp(H), Tirr(H) and Tc(H), respectively, in Figs. 1 to 3. Nelson - Le DoussaP"* line has been 
drawn following the schematic plot in Fig. 1 of Ref. 24 and made to coincide with the 
disappearance of PE. The inset shows the loci of peak effect temperatures Tp(H) at low fields 
(H < IkOe; // c) in crystals A, B and C of 2H-NbSe2. ExperimentaUy determined'^ Hci(T) 
curve for H//c has also been drawn for the sake of completeness. 
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